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ABSTRACT 

P-t 

Context, Ground-based observers have on 5-6th June 2012 the last opportunity of the century to watch the passage of Venus across the 
solar disk from Earth. Venus transits have traditionally provided unique insight into the Venus atmosphere through the refraction halo 
that appears at the planet's outer terminator near ingress/egress. Much more recently, Venus transits have attracted renewed interest 
because the technique of transits is being succesfully applied to the characterization of extrasolar planet atmospheres. 
Aims. The current work investigates theoretically the interaction of sunlight and the Venus atmosphere through the full range of tran- 
' sit phases, as observed from Earth and from a remote distance. Our model predictions quantify the relevant atmospheric phenomena, 

, thereby assisting the observers of the event in the intepretation of measurements and the extrapolation to the exoplanet case. 

Methods. Our approach relies on the numerical integration of the radiative transfer equation, and includes refraction, multiple scat- 
i-C 1 tering, atmospheric extinction and solar limb darkening, as well as an up-to-date description of the Venus atmosphere. 

Results. We produce synthetic images of the planet's terminator during ingress/egress that demonstrate the evolving shape, brightness 
and chromaticity of the halo. Our simulations reveal the impact of micrometer-sized aerosols borne in the upper haze layer of the at- 
mosphere on the halo's appearance. Guidelines are offered for the investigation of the planet's upper haze from vertically-unresolved 
photometric measurements. In this respect, the comparison with measurements from the 2004 transit appears encouraging. We also 
show integrated lightcurves of the Venus-Sun system at various phases during transit and calculate the respective Venus-Sun integrated 
transmission spectra. The comparison of the model predictions to those for a Venus-like planet free of haze and clouds (and therefore 
a closer terrestrial analogue) complements the discussion and sets the conclusions into a broader perspective. 

^ . Key words. Venus transit - refraction - halo - Earth-sized extrasolar planets 
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l/"V 1. Introduction An additional reason for pursuing planetary transits has 

emerged recently. The discovery of planets in orbit around stars 

Q\ ■ This year, from about 22:10UT June 5 to 04:50UT June 6, ob- other than our Sun has opened up a ne w and rapidly-growing 

O servers across the Earth's illuminated hemisphere have the op- field in astrophysics (Mayor & Queloz, I1995K Some of those 

CN . portunity to watch Venus passing in front of the solar disk. The planets happen to periodically intercept the Earth-to-star line of 

1 transits of Venus visible from Earth are regularly spaced in inter- sight, thereby causing a dimming in the apparent stellar bright- 

^s! I vals of 8, 105.5, 8 and 121.5 years. After the June 2004 transit, ness. In combination with some form of spectral discrimination, 

• i-h ■ this year's event is the last one before 21 17. the so-called technique of transits may provide valuable insight 

" . i i i • i ■ r into the composition and state of the planets' atmospheres (e.g. 

r , ^ ■ Venus transits occupy a notable place in the history of astron- „ „ „ . i ^mz-v i , ~ ^, . . T n i <-< 2 

, . , J r „ . „ „ , J . , Seager & Deming, 2010 and refs. therein). In our Solar System, 

rrt omy and planetary sciences. Following Halley s suggestion that . °. . f Jz . . , , _ , ., r 

. . . ■ , r J , . ., ° , t ~. , , only Mercury and Venus are observable from Earth while transit- 
Venus transits might provide a measure of the terrestrial paral- • , n n- »i . . • .1 ir 
, , .? , ,„. • i i j z. ^ ing the Sun. Since Mercury s atmosphere is tenuous, the Venus 
lax, the events of the 18th and 19th centuries helped estimate the ° . ,,„„. , f r 
„ , „ „ , , , transits of 2004 and 2012 represent unique occasions for test- 
Earth-Sun distance. Also, it is often quoted that the 1761 tran- ■• • T-1J1 . • 

., , c ., c \j ■ /t • i ing the technique of transits on an Earth-sized planet hosting a 

sit provided the first evidence for the Venus atmosphere (Link, ,° , . ,-, r ^ , ■ r 

1 r./r.i t i j i i i dense atmosphere, a maior goal in the near-luture exploration of 

1959). Lomonosov correctly attributed the halo that appears at , , J ° r 

rr , ., . . . . , ■ , . r Ut extrasolar planets. 

the planet s outer terminator during ingress/egress to sunlight r 

rays refracted towards the Earth on their passage through the The interest raised by the 2004 and 2012 Venus transits is ap- 

Venus atmosphere. The subsequent transits of 1769, 1874 and parent (Ambastha et al., 2006; Ehrenreich et al., 2012; Hedelt et 

1882 provided additional opportunities for the investigation of al.. 1201 11 Kopp et al.. 120051 Pasachoff et al.. 1201 11 Schneider et 

the phenomenon and, in turn, of the Venus atmosphere. Since al., 2004 2006; Tanga et al., 2012). An ambitious international 

then, a number of theories have been issued to interpret the halo plan of observations from both Earth-based and space-borne ob- 

and its connection with the Venus atmospheric structure. A dis- servatories is underway for the 2012 event (Pasachoff. 120121 ). It 

cussion of them and of other phenomena related to refraction in thus seems appropriate to review the theory of planetary tran- 

planetary atmospheres is presented by Link ( 119691 ). sits in its specific application to the Venus transit, taking ad- 
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vantage of recent progress in the characterization of the Venus 
atmosphere and of modern computational capacity. 

This work addresses the interaction of sunlight and the Venus 
atmosphere during transit. Descriptions of the formulation and 
the prescribed model parameters are given in $2] Section ^in- 
vestigates the shape, brightness and chromaticity of the halo. 
Section $4] establishes the connection between the halo and the 
in-transit signature of the Venus atmosphere, and presents the 
lightcurves of the Venus-Sun system as observed from Earth 
and remotely. Section §|5] summarizes the main predictions. 
Ultimately, it is hoped that our model predictions can serve as 
guidelines in the interpretation of the observations to come. 

2. Refractive theory of the Venus transit 

Lunar eclipses and planetary transits are conceptually simi- 
lar phenomena that admit a common formal treatment (Link, 
1969 ). Garcia Munoz & Palle (1201 11 1 revisited the theory of lu- 
nar eclipses. Their formulation was subsequently used in the 
interpretation of a partial lunar eclipse (Garcia Munoz et al., 
1201 11 ). and to investigate the impact of refraction on the in-transit 
signature of Earth-like extrasolar planets (Garcia Munoz et al., 
120 12i We now apply the formulation to the Venus transit. Since 
the refractive theory of planetary transits has thus far received 
marginal attention, we expend some effort describing the model 
approach. 



2.1. The Radiative Transfer Equation (RTE) 

The steady-state RTE for a refractive medium in the geometrical 
optics limit has been presented by, e.g., Zheleznyakov (1967) 
and Enome ( 1969 ): 



s- V 



' L(x,s) 
v n 2 (x) 



L(x,s) 7(x,s) 
n £ {x) n £ (x) 



(1) 



Here, L(x, s) is the radiance at x in direction s, y(x) and n(x) 
stand for the extinction coefficient and index of refraction of the 
medium, respectively, and J(x, s) is the source term, which may 
include separate terms for scattered radiation and emission from 
within the medium. L(x, s), y(x), n(x) and other associated mag- 
nitudes are wavelength dependent, which should be explicitly 
stated in the notation. For simplicity, however, that characteris- 
tic is omitted throughout the text. The boundary conditions for 
the direct and diffuse components of Eq. (Q]) are introduced in 
j32.2l and [2.3l respectively. 

For n(x) constant, Eq. (Q]( turns into the well-known RTE 
for a non-refractive medium. In that case, the photons follow 
straightline trajectories between scattering events. For a non- 
absorbing, non-emitting, non-scattering medium (y(x)=J(x)=Q) 
with possibly spatially-dependent «(x), Eq. ([1) reduces to the 
conservation of L(x, s)/n 2 (x) on the refraction-bent ray trajec- 
tories. That magnitude is sometimes referred to as the Clausius 
invariant, and generalizes the conservation of radiance in media 
of constant index of refraction. 

As in radiometry of non-refractive media, the irradiance, or 
net flux of light that crosses an elementary surface of normal 
vector n x<) (not to be mistaken for the index of refraction at 
xo, «(xo)) from directions So, at the observer's location xo is 
obtained by integrating the directional radiance L(xo, So) over 
solid angle: 



F(*o) 



L(x ,s )s ■ n Xo dQ(s ). 



(2) 



In the general application of Eq. (f2]i to the Venus transit, the 
integration domain must contain the solid angles subtended from 
xo by planet and atmosphere together, dClo, and the Sun, <9Q . 

For the numerical evaluation of Eq. (f2|, it is convenient to 
work with auxiliary variables. Concentrating for now on the in- 
tegral over dQcp, introducing: 

/j. = s ■ n x<) , 

and the azimuthal angle <f>, using: 

c/f2(So) = dfidip, 

and the normalized variables: 

<? = (m 2 - - a4) and n = 0/*. 

Eq. (f2j) can be rewritten as: 

F(xo)=7r(/4~/4) f f mrfidtdr}, 
Jo Jo 

where: 



(3) 



i4a - 1 - (fmA/do) 2 , 
H 2 m = 1 - (R 9 /d Q ) 2 , 
rroA = + hjoA, 

is the Venus mean radius, /ztoa is the top of the atmosphere 
altitude, and do is the observer-to- Venus centre distance. We 
may also define the impact radius r b for each fi 2 through: 

fx 2 = 1 - (r b /do) 2 . 

Here, r/, has the meaning of the distance of closest approach to 
the Venus centre for straightline rays traced sunwards from xo- 
Essentially, Eq. (0 maps the original integral onto a plane of 
incident directions So={£, rf\ sampling the field of view from xo- 
Figure 1 in Garcia Munoz & Palle ( 1201 II ) sketches the meaning 
of some of the variables introduced above. 

2.2. The direct component of sunlight 

The usual treatment of the RTE separates the direct and diffuse 
components. Neglecting the emission from within the medium 
and taking J(x, s)=0, one obtains for the direct radiance: 



L D (*o, So) L D (Xq, Sq) 



n 2 (x Q ) 



n 2 (x G ) 



'"I 



exp[- I y(s)ds], 



(4) 



for ray trajectories F connecting xo and the solar disk. Here, s 
is the path length on the ray trajectory and x and s Q are the lo- 
cation and direction of departure of the sunlight ray on the solar 
disk. The squared refraction indices at the observer's site and at 
the Sun's photosphere, n 2 (xo) and n 2 (x ) respectively, can be 
taken as a; 1 in Eq. (0]i at visible and near-infrared wavelengths. 
Ld(x q , s ) is the radiance emitted by the Sun and enters the for- 
mulation as a boundary condition. We assume that it depends on 
the emission angle, cos~'(/L(oX through a polynomial of the form 
U= Xli=o M 'A'o' an d tnat tne wavelength-dependent coefficients 
Ui account for the chromatic effects in the solar output. A good 
representation of the solar output near the Sun's edge, ju — >0, is 
important for the simulation of the ingress/egress phases of the 
transit. The adopted wavelength-dependent limb-darkening co- 
efficients Ui are from Pierce & Slaughter ( 1 19771 ) and Pierce et 
al. (119771) . Generally, the coefficients are dependent on the star's 
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metallicity, a fact that must properly be accounted for in the anal- 
ysis of extrasolar planet lightcurves (Claret, 20001 Sing, [2010). 
Thus, Lo(x Q , s q )=B q U, where B G is the stellar radiance at the 
Sun's center, // Q =l. (Note also that Y^=o «i=L) Figure ([TJ shows 
the implemented t/(yU©) at selected wavelengths. 
Considering the above, to a good approximation: 

L D (x , s ) = B Q U exp [- j ' y(s)ds], (5) 

which is the Beer-Lambert law of extinction in a non-refractive 
medium. Using Eq. (0 does not mean that refraction plays a mi- 
nor role in the formulation. Refraction enters Lo(xo, So) at the 
time of devising the ray trajectory T. For that purpose, we use 
the tracing scheme based on path length by van der Werf (2008) 
and trace each ray trajectory from xo sunwards after specify- 
ing the Venus atmospheric profile for the wavelength-dependent 
refraction index. 

For the evaluation of Eq. (O over 3f2cp, the ^-77 field of view 
at xo is subdivided into NgXN, ; (=200x360) evenly-sized bins, 
each of them associated with an incident direction So- For each 
direction So={£, rf\ a ray trajectory F is traced, and the corre- 
sponding optical opacity going into Eq. (f5J) is calculated. The 
rays failing to connect xo with the solar disk have Lb(xo, So)=0. 
In the halo, during ingress/egress, Lb(xo, So)=0 occurs for rays 
crossing nearly unrefracted through the upper layers of the at- 
mosphere, but also for rays passing through the lowermost al- 
titudes and whose trajectories become deflected away from the 
Sun. During transit, the trajectories failing to connect xo with 
the solar disk impose a ring of altitudes at the planet's termina- 
tor that cannot be probed from a remote distance. The latter is 
discussed in detail by Garcia Munoz et al. (2012) and will be 
noted again for the specific case of Venus in 33.11 and 14.11 By 
going through the entire N^xN l; matrix of directions So, the im- 
age of the solar disk at xo emerging from the planet's terminator 
is constructed. 

We could adopt the same strategy to evaluate Eq. ([3]) over the 
rest of dQ. not contained in dQc^. This is, however, rather ineffi- 
cient given the large ratio of Sun/Venus radii, R e /Rq. Instead, 
we integrate analytically the irradiance over the unobstructed 
view of the star, FQ=7i(R Q /l Q ) 2 Yj] = Q2ui/(i + 2), where l Q is the 
observer-to-Sun centre distance, and subtract the part corre- 
sponding to d£lq calculated by Monte Carlo integration. 

Taking the observer's location xo for carrying out the inte- 
gration of Eq. (f3j) for the direct sunlight component is in line 
with the early formulations of the photometric theory of lunar 
eclipses (Link, IT9621 [1969b . Link ( [19691 ), in Section 1.2.21 of 
that monograph, comments on the equivalence of that approach 
with other formulations that alternatively carry out the integra- 
tion at the solar plane. Link (1969) also refers to the so-called 
attenuation by refraction that occurs by the compression of the 
solar disk image at the planet's terminator during ingress/egress, 
and notes that its effect is implicitly accounted for in the projec- 
tion of c/£2(so) on the solar disk. Attenuation by refraction occurs 
always that sunlight rays traverse a planet's stratified atmosphere 
(Hays & Roble, |T9681 ), and merely represents the change in cross 
section undergone by a pencil of sunlight rays upon crossing the 
stratified medium. 

Figure (0 offers some visual insight into this point. The rays 
in the sketch are evenly spaced in angle at their departure from 
xo towards the Sun. The rays gradually separate as they tra- 
verse the atmosphere due to differential refraction between al- 
titude layers. Assuming that the atmosphere is transparent and 
solar limb darkening is omitted, Eq. ([3]) would simply lead to 



a measure of the halo's angular size perceived at xo- Since in 
such conditions the radiance is very approximately conserved 
on ray trajectories, the contribution of a solar disk surface ele- 
ment is ultimately limited by the size of its refracted image at the 
planet's terminator. Sidis & Sari (2010) set out from this conclu- 
sion to evaluate the halo's overall brightness for transiting extra- 
solar giant planets with transparent atmospheres. The factor by 
which the pencil's cross section is compressed depends on both 
the Sun- Venus and Venus-observer distances as well as on the 
atmosphere's refractive properties. Figure 0} sketches the inte- 
gration of the direct radiance over a differential element of solid 
angle at the observer's site. 

2.3. The diffuse component of sunlight 

A fraction of the sunlight photons incident on Venus is redirected 
in the forward direction after one or more scattering collisions 
with the medium. To estimate the contribution of diffuse sun- 
light, we must solve the RTE including: 

7(x,s) =/?(x) f p(x,s,s')£(x,sVQ(s') (6) 

as a source term of the equation, with B{x) being the medium's 
scattering coefficient and p(x, s, s') the normalized phase func- 
tion for changes s'— >s in the direction of light propagation. The 
given 7(x, s) is identical to the form of the scattering term in 
non-refractive media (Ben-Abdallah et al., 2001; Kowalski & 
Saumon, 2004; Fumeron et al., [2005 ). The integration over the 
part of dQ. not included in dQ.^ is unnecessary if the top of the 
atmosphere is placed high enough. 

Working in a refractive medium presents two additional sub- 
tleties in the treatment of the diffuse RTE. First, the radiance is 
affected by the reciprocal of n 2 (x), although we are allowed to 
take n(x)»T at the relevant wavelengths. Second, the photon tra- 
jectories in between scattering events follow refraction-bent tra- 
jectories. However, scattering itself alters the direction of light 
propagation by angles likely to exceed the deflection caused by 
refraction. We may therefore also omit this aspect, and treat the 
diffuse RTE as in a non-refractive medium. 

The practical solution to Eqs. (Q} and (O is carried out by 
Mont e Carlo integration (O'Brien, [19921 Garcia Munoz & Palle, 
1201 It with the proper expression for F G at the Venus orbital dis- 
tance as a boundary condition. The scheme is nested in the si- 
multaneous integration of Eq. ([3]) over the Venus disk. 

2.4. The reference Venus atmosphere 

Traditionally, observations of Venus during ingress/egress have 
provided valuable information on the tilt of the planet's rota- 
tional axis and equator-to-pole variations in the atmospheric 
structure (Link. fT959b . Thus, for the simulations we require a de- 
scription of the atmospheric optical properties on a global scale 
as well as the latitudinal dependence for some of them. 

For the temperature, we adopted the Seiff et al. ( 1985) pro- 
file for a latitude of 45° up to 100 km and the Hedin et al. ( 1983 ) 
noon profile there upwards. The atmosphere was assumed to be 
composed of CO2-N2 in a 96.5-3.5% proportion, and the back- 
ground density was obtained by integrating the hydrostatic bal- 
ance equation. 

At the altitudes at which the atmosphere becomes optically 
thick, the opacity is largely dictated by aerosols rather than by 
the atmospheric gas. Thus, latitudinal variations in the gas profile 
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are expected to be of less importance than the respective varia- 
tions in the aerosol loading. Similarly, the errors introduced in 
the refracted angles from neglecting the latitudinal variation in 
the density profiles are minor. 

The gas refractivity is needed to trace the ray trajectories and 
to calculate the Rayleigh scattering cross sections. We adopted 
for C0 2 the fit between 0.48 and 1.82 jum by Old et al. JTWTT) , 
and for N 2 the fits between 0.2 and 2.1 fim by Bates ( 11984| >. 
The fits were truncated at 2 jjm at the longer wavelengths and 
extrapolated at the shorter wavelengths. 

The Venus atmosphere is heavily loaded with aerosols that 
become optically thick in the nadir at about 70 km. By con- 
vention, aerosols above and below the optical thickness T na dir=l 
are usually referred to as haze and cloud particles, respectively 
(Esposito et al., I19831 I. Limb opacities are enhanced by a geo- 
metrical factor (Barth, [19691 of ~ ^2ttR^/ H da ~l00, estimated 
with an aerosol scale height H- da -~4 km, which means that we 
will hereon only consider the aerosol opacity due to haze. 

To calculate the haze opacity we followed the prescriptions 
by Molaverdikhani et al. (2012), which are consistent with those 
by Pollack et al. (119801 1. Thus, from 60 to 80 km, we adopted log- 
normal size distributions of aerosols that depend on the particle's 
radius r a through: 

11 , , 

n a (r a ) = — — exp (-(In r a - In r g ) /2 In cr g ). 

V27iTn cr g r a 

Mode 1 haze is characterized by r g =0.15 /urn and <x g =1.5 (ef- 
fective radius and variance of r e ff=0.23 fim and v e ff=0.18, re- 
spectively), and a number density profile that drops with the 
background pressure as 900xp/p(60 km) cm 4 . For mode 2 
haze, r g =1.0 /mi, cr g =1.21 (r e ff=1.09 jum and v e ff=0.037), and 
the number density profile is 120xp/p(60 km) cirT 3 . For our 
simulations, we extended indefinitely upwards the p/p(6Q km) 
profile. Indeed, mode 1 particles are long known to exist above 
80 km, and two works (Wilquet et al., |20U9l de Kok et al.. [20TTb 
recently identified sizeable amounts of mode 2 particles up to 
90-95 km. Both r e jf and v e ff are defined by Hansen & Travis 
d!974| l. Note, though, that the definition of <r g in that work corre- 
sponds to In cr g here. We assessed that the effect of extending the 
aerosol profile above 100 km was negligible because the aerosol 
layer becomes effectively thin at such altitudes. The wavelength- 
dependent cross sections, albedos and scattering phase functions 
for each mode are calculated from Mie theory (Mishchenko et 
al., [2002 ). The complex refractive indices for the aerosol parti- 
cles (H2SO4/H2O solutions at 84.5% by weight, Molaverdikhani 
et al. (120121) ) were borrowed from Palmer & Williams d 1 975b . 
Figure (O graphs the extinction profile at 60 km from 0.2 to 5 
yt/m, that appears to be dominated by mode 2 particles. Figure 
(0) shows the scattering phase function of each mode at specific 
wavelengths. 

The aerosol profile based on Molaverdikhani et al. (120121 ) 
reaches T na di r =l at 1.6 pm near 70 km, which may represent 
a globally-averaged cloud top altitude. A limb opacity Tu m b=l 
is reached about 4-5 scale heights higher and therefore likely 
above 85 km. This globally-averaged description is correct to 
first order, but the appropriate simulation of the halo requires 
a more accurate characterization of latitudinal variations in the 
cloud top level. Recent altimetry of the Venus clouds from Venus 
Express shows that T na di r = 1 is reached at about 74 km near the 
equator but may drop to 64 km at the poles, with transition zones 
on both hemispheres at about 60° latitude (Ignatiev et al., 2009). 
To account for this latitudinal dependence, we shifted at each 
latitude the aerosol profiles upwards or downwards as needed to 



match the cloud top altitudes graphed in Fig. 8b of Ignatiev et al. 
(I20091) . 

3. The halo 

The halo (or aureole) is the refracted image of the Sun that forms 
at the planet's outer terminator for small Sun- Venus angular sep- 
arations, with possibly a minor contribution from diffuse sun- 
light. Its appearance is determined by the geometrical configu- 
ration of the two objects relative to the observer and by the op- 
tical properties of the Venus stratified atmosphere. This section 
addresses the visual phenomenon as contemplated from a geo- 
centric vantage point for the conditions of the June 2012 transit. 
The geometrical parameters needed in the simulations are either 
adopted or inferred from the JPL HORIZONS service (Giorgini 
et al.. 1 19961 ) or the New Horizons Geo Viz software Q 

3.1. Some geometrical considerations 

Figure (|6]) shows a sequence of simulations of both the ingress 
and the egress. Each simulation corresponds to a value for the / 
factor introduced by Link ( 1969 ) and used by Tangaet al. ( 120121 ). 
This factor stands for the fraction of the planet's diameter along 
the Venus-Sun line of centres beyond the Sun's limb. Thus, /= 
1, 0.5 and are appropriate for external contact, mid-ingress/- 
egress and internal contact, respectively. For a better visualiza- 
tion, the atmospheric ring from 50 to 100 km altitude is stretched 
by a fixed factor in all images. Ignoring for the time being the 
sunlight extinction caused by the gas and haze, the Sun's image 
at the planet's outer terminator appears notably diminished with 
respect to the apparent size of the unocculted solar disk. The 
diminishment constitutes the visual representation of the attenu- 
ation by refraction. 

Referring to Table ([1), the deflected angle for rays cross- 
ing the atmosphere with an impact distance ri,-Rq=55 km is 
ffrefr^l , thus exceeding the solar angular diameter as viewed 
from Venus (~2/? G /0.726 AU~0.734°). This means that, in ge- 
ometrical terms, the halo representations of Fig. (O are the re- 
fracted images of the entire solar hemisphere visible from Venus. 
The innermost/outermost contours are the projections of the far- 
ther/closer edges of the solar disk, respectively. 

The bottom left and top right images of Fig. (O show the 
planet shortly after contact II and before contact III, respectively. 
At those instants the atmosphere does not form a refracted im- 
age at the inner terminator below ~76 km on the axis between 
the centres of the two objects. At the outer terminator, a refracted 
image occurs only above ~59 km. This asymmetry is caused by 
the unequal distance from each terminator to the opposite solar 
limb. The oval ring within which a solar image is not formed dic- 
tates the range of atmospheric altitudes that cannot be accessed 
from Earth regardless of atmospheric opacity. This refraction ex- 
clusion ring becomes circular as the planet moves towards mid- 
transit. 

3.2. Shape and brightness of the halo 

The simulations of Fig. (O assume a continuum wavelength of 
0.55 pm and the latitude-dependent extinction coefficient de- 
scribed in 32.41 To take into account the varying altitude of 
clouds along the terminator, we solved the RTE over the entire 
Venus disk for a range of /z c i ou d values and, subsequently, com- 



H. Throop, http://soc.boulder.swri.edu/nhgv/ 
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posed the displayed images by assembling the relevant latitude- 
^cioud sectors. 

Venus enters the transit with its rotational axis tilted by about 
45° clockwise with respect to the Venus-Sun line of centres. At 
egress, the tilt is about 60° and anticlockwise. The tilts were esti- 
mated graphically from Geo Viz simulations. In combination, the 
inclination of the rotational axis and the latitude-dependent alti- 
tude of the cloud tops break the symmetry of the halo expected 
from purely geometrical arguments. The asymmetry is clearly 
discerned in the contours of Fig. ©, each of them correspond- 
ing to a value for <!>=- log 10 (L D (x o , s )/B e ). 

At ingress, the earliest evidence of the halo occurs over the 
Venus North Pole and appears detached from the solar disk until 
about mid-ingress. The bright cap finally turns into a full ring 
that encircles the planet before contact II. Essentially, the egress 
reproduces the ingress sequence in the reverse order. During 
egress, the larger tilt of the rotational axis causes the halo to 
appear attached to the solar disk for a longer duration. 

Figure © provides insight into the instantaneous vertical ex- 
tent of the halo. Importantly, the halo's outermost contour re- 
mains confined to altitudes less than 100 km for most of the 
ingress/egress. The statement is easily verified with a quick es- 
timate specific to the instant halfway through ingress/egress. At 
that moment, a ray emanating from the nearer solar limb and 
propagating on the azimuthal plane that includes the Venus- 
Sun line of centres must be deflected by an angle ~R^/0.126 
AU+Ry/0.214 AU~42 arc sec to reach the geocentric observer. 
Refraction yields such a deflection at ~80 km altitude, as seen in 
Table ([1]). Thus, at mid-ingress or -egress only rays traced from 
xo on the plane containing the line of centres with less than 
that may connect the observer and the solar disk. 

The conclusion serves to highlight the importance of haze on 
the halo's changing appearance because the atmosphere is likely 
optically thick in limb viewing at those altitudes. The halo's 
inner contour is largely dictated by extinction in the optically 
thicker layers of the atmosphere. Overall, the altitudes contribut- 
ing effectively to the halo range over 1-2 scale heights. 

3.3. The imaged solar disk 

It is interesting to trace the rays forming the halo back to the 
solar surface elements they departed from. In geometrical terms, 
the halo is the refracted image of most of the solar disk visible 
from Venus or Earth. However, haze extinction concentrates the 
rays effectively contributing to the halo's brightness onto a nar- 
row strip near the solar limb. Our ray tracing scheme does indeed 
show that during ingress/egress the halo is essentially formed 
by the solar disk occulted behind Venus. That solar region is 
severely limb-darkened, which adds a significant radial gradient 
to the solar output brightness. In the serendipitous situation that 
a dark sunspot is occulted by the planet during ingress/egress, 
one may expect that the sunspot will turn up as a dimmer halo. 

3.4. Chromaticity of the halo and limb-integrated lightcurves 

The halo's shape and brightness depend on the atmospheric re- 
fractivity, sunlight extinction and solar limb-darkening. In turn, 
these factors depend on wavelength. Focusing on the visible and 
near-infrared regions of the spectrum, we investigated them at 
0.43, 0.55 and 0.8 /mi, which roughly match the center wave- 
lengths of the B, V and I filters, respectively. 

The refractivity of the Venus atmosphere varies with wave- 
length, which may introduce a relative displacement between the 



images formed at different colors. Considering again Fig. ©, we 
can devise two rays of wavelengths Ai and A2 departing from a 
common solar surface element, both of them reaching the ob- 
server. Since the overall deflections are nearly identical, we can 
readily estimate their altitudes of closest approach in the atmo- 
sphere. According to the approximate, analytical formula given 
by Baum & Code ( 1953) for the refraction angle in an isothermal 
atmosphere, the number densities at closest approach of the two 
rays are in the ratio ni/«2~v '\J\ % , with v being the refractivity at 
the specified wavelength. For Venus, this means that if ^1 =0.43 
and /l2=0.8 fim the rays have their closest approach at densities 
in the ratio «i/«2~0.9747. For an atmospheric scale height of 4 
km, this means that the altitudes of closest approach differ by 
about one hundred meters. Thus, the displacement of the images 
due to refraction is relatively minor. 

Haze dominates the atmospheric opacity from the UV to the 
NIR at the altitudes of formation of the halo. Recalling from Fig. 
(O the spectral dependence of y(x), it is apparent that the sun- 
light emerging from the Venus terminator must exhibit a mod- 
erate red tone imprinted by haze extinction. The reddening is, 
however, much less than expected for a Rayleigh atmosphere, 
in which case y(x)~/T 4 . Sunlight emitted from the solar limbs 
is eminently red, as seen in Fig. (Q]i. Thus, solar limb darkening 
must also contribute towards the reddening of the halo. 

To study this in a quantitative manner, we produced radially- 
averaged radiances of the halo as follows: 



^ L d (Xq, Sq) ^ _ 



jL D (x ,So)/B Q d(M 2 /2) 1 
J d(fi 2 /2) 



d? (7) 



which is equivalent to averaging over altitude from the ground 
to 200 km or over 1 arcsec for ground-based observations. Thus, 
<Lo(xo, So)/B©> is a valid measure of the vertically-unresolved 
brightness with respect to the brightness at the Sun's centre. 

Curves of the radially-averaged radiance for A-0.43, 0.55 
and 0.8 /mi during ingress and egress are graphed against the 
azimuthal angle (f> in Fig. (0. The curves show a strong time- 
dependence and reveal with an increased brightness the lower 
clouds occurring at the Venus North Pole. The color dependence 
of the lightcurves near the Sun- Venus line of centres is due to 
both solar limb darkening and atmospheric sunlight extinction. 
Close to the solar disk edge (the left and right pedestals of the 
figure), the higher atmospheric altitudes contribute significantly 
and only solar limb darkening introduces a color dependence. 

These curves constitute a direct target for comparison with 
photometric measurements. Indeed, Tanga et al. (2012) recently 
published limb-integrated curves from observations of the 2004 
transit. Figure ([8]l shows with symbols their G-band lightcurves 
obtained during egress at the Dutch Open Telescope on La Palma 
as digitized from their Fig. 8. Relative intensities are represented 
against the azimuthal angle and reveal the expected U shape as 
well as a peak in brightness displaced from the pole. The solid 
curves are our model lightcurves at 0.43 /mi for the 2004 egress 
conditions. Tanga et al. (2012) normalize their lightcurves with 
the solar disk brightness at about one Venus radius from the solar 
edge. For consistency, we re-scaled our radially-averaged radi- 
ances by the appropriate limb-darkening factor, so that the two 
sets of curves become directly comparable. The comparison for 
three different instants during egress is overall satisfactory, espe- 
cially taking into account that we assumed a global description 
for the cloud top altitude. Importantly, to achieve such an agree- 
ment, we needed to offset the location of the minimum cloud top 
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altitude by about 20° away from the pole. If confirmed, that fea- 
ture might indicate that occasionally the cloud top altitude does 
not descrease monotonically with latitude towards the poles. 

Following the comparison with Tanga et al. (2012), that work 
reports that as factor / increases the halo brightness decays more 
rapidly at red wavelengths than at visible and ultraviolet wave- 
lengths. As discussed earlier, our model predicts an overall red 
halo, which is somewhat at odds with the Tanga et al. (2012) 
findings. At any rate, the amplitude of the error bars in their ex- 
periment might also accomodate a small reddening. New obser- 
vations of the 2012 transit will hopefully provide a better insight 
into the issue. 

One may consider to use model lightcurves such as those 
in Fig. (jTj to fit the empirical lightcurves and retrieve cloud al- 
titudes along the planet's terminator. Tanga et al. (120121 1 have 
already approached the task. To understand to what extent the 
model lightcurves depend on the prescribed atmospheric opti- 
cal properties, we produced limb-integrated lightcurves with the 
cloud top altitude shifted at all latitudes by + 1 km with respect 
to the reference cloud altimetry. The resulting lightcurves are 
given in Fig. (0 and show that shifting the cloud tops by ±1 km 
leads to variations in the limb-averaged radiances by factors of 
~2 or more. Ultimately, the precision of the method to retrieve 
the cloud top altitude will depend strongly on the quality of the 
photometric data. Spatial resolution in the data is not anticipated 
to be critical, since resolutions of about 1 arcsec are commonly 
achieved also from ground-based observatories. A cadence of 
one image per minute or less suffices to reveal the temporal evo- 
lution of the lightcurves. 

Some observers of the 2012 Venus transit may also attempt 
the spectroscopic characterization of the halo. The likely result 
of such an experiment will be similar to a transmission spectrum 
of the Venus atmosphere from an altitude of 85 km or higher 
with a continuum dependent on wavelength through both haze 
extinction and solar limb darkening. 

3.5. Diffuse sunlight contribution to the halo 

One more issue to explore is the contribution of diffuse sun- 
light to the overall brightness of the halo. In between external 
and internal contacts, the Sun- Venus-Earth angle is more than 
~179°, which favors forward scattering from the planet's termi- 
nator. Most atmospheric aerosols scatter efficiently in the for- 
ward direction, with the scattering efficiency depending strongly 
on the ratio x e ff=2nr e ff/A of the effective radius for the particle 
size distribution and the incident sunlight wavelength (Hansen 
& Travis, 1974). As seen in Fig. (0, the peak in the scattering 
phase function becomes enhanced at the shorter wavelengths and 
is particularly stronger for the larger mode-2 haze. 

Figure ( fTOb graphs the multiple-scattering diffuse radiance 
at mid-transit (with the Sun, Venus and Earth perfectly aligned) 
against the impact radius as calculated with the Monte Carlo 
method of 32.31 For simplicity, we assumed that the cloud top 
lies at a constant altitude of 70 km at all latitudes. The radiances 
peak at ~88 km, that corresponds approximately to the rii m b=l 
level. Shifting upwards/downwards the clouds by a few kilome- 
ters merely shifts the radiances by the same distance but has a 
negligible impact on the vertically-integrated magnitudes. 

Integration of the curves in Fig. ( fTOb results in radially- 
averaged diffuse radiances <L ( ;(xo, So)/fi©> of ~3xl0~ 5 , 
2xl0~ 5 and lxl0~ 5 at 0.43, 0.55 and 0.8 yum, respectively. 
Comparison with the lightcurves of Fig. (|7]) indicates that re- 
fracted sunlight dominates the halo's brightness over most of the 
ingress/egress. 



4. Venus-Sun disk-integrated lightcurves 

We have thus far assumed that the Venus disk is spatially resolv- 
able from Earth, a necessary condition to investigate the connec- 
tion between the halo's structure and the cloud altimetry along 
the planet's terminator. However, the question of whether the at- 
mospheric properties of a transiting planet can be disentangled 
from the planet-star blended signal has become key in the re- 
search of extrasolar planets (Seager & Sasselov, 120001 Brown, 
I200T1 Hubbard et al.. l200TT) . 

The so-called technique of transits aims to identify the sig- 
natures of the planet and its atmosphere from comparative mea- 
surements of the planet-star light out of and in transit. In-transit 
spectroscopy (or spectrophotometry) is providing valuable in- 
sight into exoplanetary atmospheres, especially for the giant 
planets that constitute most of the exoplanet detections to date. 
As the sizes of discovered extrasolar planets shrink and technol- 
ogy improves, the prospects for characterizing Earth-sized plan- 
ets with this technique, including Venus-like ones, look closer. 
Next, we describe the Venus transit as if observed in integrated 
light from a terrestrial distance and from a remote vantage point. 



4.1. Disk-integrated lightcurves from Earth and far away 

Figure (fTTT) shows model transit depths at 0.55 /im for the Venus- 
Sun system as observed from Earth. The transit depths are ex- 
pressed in units of irradiance for the unobstructed Sun. The 
overall shape of the lightcurve is determined by the ratio of 
Venus/Sun solid angles, ~l/30, and the limb darkening function 
at the specified wavelength. 

Refraction introduces unique features in the lightcurves near 
external and internal contacts. In Fig. (JTTJ, we show also the 
difference in transit depths when refraction is taken into ac- 
count and when it is omitted. If the Venus atmosphere was free 
of haze and cloud particles, the differential transit depth curve 
would show that the planet's size is determined by a combina- 
tion of Rayleigh scattering and refraction. A refraction exclu- 
sion ring appears at altitudes from ~59 to ~76 km, depending 
on the distance from the planet's local terminator to the solar 
limb, and contributes most significantly to the disk-integrated 
lightcurve near internal contact. Note also that the depth varies 
by ~3xl0~ 7 from mid-transit to internal contact, which means 
that the planet's measurable radius apparently expands by a few 
km. Out of transit, the transit depth inverts its sign because sun- 
light refracted by the planet contributes positively to the bright- 
ness of the Venus-Sun system. The magnitude of this additional 
brightening amounts to ~7xl0~ 7 at external contact. Including 
haze in the calculation of the optical opacity reduces by orders 
of magnitude the effects of refraction in the differential transit 
depths, as the dotted curve in the figure indicates. At any rate, 
such small perturbations fall well below typical irradiance fluc- 
tuations due to solar oscillations over periods of minutes (Kopp 
etal.. l2lX)5l 

The disk-integrated lightcurves measured by a remote ob- 
server would exhibit the same phenomena described above. The 
considerably smaller ratio of Venus/Sun solid angles in that case, 
~Rq/R ~l/115, would cause the same phenomena to appear 
about one order of magnitude fainter in the combined planet-star 
signal. Figure ( fT2l demonstrates the disk-integrated lightcurves 
in the limit of an infinitely distant observer. 
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4.2. Disk-integrated transmission spectra 

The monochromatic lightcurves presented in ^4.11 might also 
be calculated at wavelengths exhibiting discrete molecular ab- 
sorption features. We produced spectra for the transit depth at 
wavelengths from about 0.3 to 5 fim. In addition to CO2, the 
atmosphere was assumed to contain small amounts of N2, CO 
and O (Vandaele et a l., [20081 Hedin et al„ [1983b and H 2 (1 
ppm; Fedorova et al.. l2008l l. The line parameters for molecular 
transitions were borrowed from HITRAN 2008 (Rothman et al., 
2009 ). The transit depths were degraded to a resolving power of 
about 10 3 and expressed in terms of the equivalent height. (The 
equivalent height is the height of an opaque slab that would pro- 
duce an identical transit depth.) In equivalent heights, the spectra 
measured from Earth and from a remote distance are nearly iden- 
tical, yet the associated stellar dimming may differ by an order of 
magnitude. Whether observed from Earth (do/l o ~0.274) or from 
a remote distance (do/le~l), the solar irradiance blocked by the 
Venus disk and atmosphere is Fd +cl ~n((Rq + h eq )/do) 2 B U([i Q ), 
where we tacitly assume that the planet is nearly point-like and 
the background solar radiance is determined by a single fi Q de- 
pendent on the planet's phase. In relative terms, the solar dim- 
ming during transit is Fd+ a IF Q . 

The curves in Fig. ( TT3b correspond to mid-transit (/z Q =l) 
spectra calculated for an atmosphere free of haze and cloud par- 
ticles and for a hazy atmosphere with cloud tops at 70 km. The 
latter seems to reproduce most of the structure of the in-transit 
spectrum presented by Ehrenreich et al. (1201 21 ) . The most ob- 
vious difference between the Rayleigh and hazy atmospheres 
refers to the continuum level, that lies considerably higher in 
the latter. For the Rayleigh atmosphere, the continuum level is 
largely determined by refraction at wavelengths longer than ~1 
yum, with Rayleigh scattering contributing shortwards of that. 
For the hazy atmosphere, the continuum is dictated by the haze 
optical properties throughout the spectrum. The simulations for 
phases other than mid-transit (not shown) indicate that for the 
Rayleigh atmosphere the transit depth of molecular bands is 
highest at mid-transit and a few kilometers less when the planet 
is near contacts. This modulation is attributable to refraction, and 
becomes entirely negligible for the hazy atmosphere. 

The comparison is relevant because to some extent a clear- 
atmosphere Venus would be akin to our own Earth. Regarding 
the molecular structure of the spectra, it is overwhelmingly dom- 
inated by rovibrational bands of CO2 . Benneke & Seager ( 20 1 2 ) 
have investigated what can unambiguously be inferred from in- 
transit spectroscopy of a super-Earth, and noted that the lacking 
evidence of Rayleigh scattering may impede the identification of 
the main atmospheric constituent. Interestingly, it appears that 
the in-transit spectrum of Venus would reveal little on the iden- 
tity of the main atmospheric constituent. 

A final comment on the SO2 molecule. The SO2 molecule 
plays a key role in the Venus photochemistry and in the forma- 
tion of H2SO4 in the clouds (Mills & Allen, [20071 . The S0 2 
molecule absorbs at wavelengths below ~0.32 //m, at the limit 
of accessibility for ground-based observations imposed by the 
terrestrial ozone cut-off at ~0.3 //m. Belyaev et al. (2012) have 
retrieved the SO2 mixing ratio for altitudes below 105 km. The 
precise amounts appear to be sensitive to the retrieval approach, 
which leads to values somewhere between 10~ 7 and 10~ 6 at alti- 
tudes above the Ti; m b=l level of ~90 km. Figure (TBI shows the 
equivalent height near 0.3 fim at high resolving power for SO2 
mixing ratios of 10~ 7 and 5xl0~ 7 . The main conclusion to ob- 
tain from the graph is that SO2 introduces a distinct signature 



that might be captured in spectroscopic observations of the tran- 
sit if enough signal-to-noise ratio is achieved. 

5. Summary 

The current work presents theoretical insight into the signature 
imprinted on sunlight by the Venus atmosphere during the pas- 
sage of the planet in front of the solar disk. The conclusions 
are specific to the June 2012 event but are easily generalizable 
to other transits. Our model approach includes refraction, multi- 
ple scattering and extinction by gases and aerosols in the atmo- 
sphere, as well as solar limb darkening, and handles all phases 
of the transit (including the out-of-transit, ingress/egress and in- 
transit stages) in a consistent manner. To the best of our knowl- 
edge this is the first theoretical investigation of the Venus transit 
at this level of detail. 

We investigated the halo that forms at the planet's outer 
terminator near ingress/egress, predicting its shape, brightness 
and chromaticity. Both the inner and outer contours of the halo 
are considerably affected by haze in the atmosphere above the 
clouds and below 100 km. The impact of haze opacity on the 
halo structure explains the variability of halo patterns from 
ingress to egress and from event to event. We estimated that scat- 
tered sunlight contributes less than refracted sunlight for most of 
the ingress/egress between internal and external contacts. Since 
our predictions are based on a realistic description of the Venus 
atmosphere it should be possible to verify them provided that ac- 
curate measurements become available. Further, our formulation 
opens a way for the characterization of the Venus atmosphere 
from observations during the transit. In retrospective, it is fair to 
state that if the Venus haze extended for 2-3 scale heights higher 
than it does, it is likely that the identification of a Venus atmo- 
sphere in the 18th century would have to have waited longer. 

We presented lightcurves for the combined Venus-Sun sig- 
nal during transit and discussed some features of the lightcurves 
associated with refraction. We also presented model in-transit 
spectra of Venus and of a haze- and cloud-free Venus-like planet. 
Haze appears recurrently as a key element in the interaction of 
sunlight and the Venus atmosphere. The disparate vertical ex- 
tent of aerosols in their atmospheres establishes a considerable 
difference between Venus and Earth. Had the proper technology 
existed, it might have been appropriate to investigate the historic 
record of observed Venus transits. In this respect, the 2004 and 
2012 events may set a valuable start point for comparisons be- 
tween transits. 

Finally, the Venus transit of 2012 arrives at a moment of great 
activity in the exploration of extrasolar planets by the method of 
transits. We have seen that high-altitude aerosols introduce spe- 
cific challenges to the characterization of the atmospheric com- 
position and cloud top altitude for a transiting Venus. Since it is 
still unknown how frequently high-altitude clouds may occur in 
the atmospheres of extrasolar planets, there are reasons to think 
that the challenges posed by Venus might also be generally en- 
countered in the characterization of those planets. Indeed, recent 
work on hot Jupiter HD 189733b (Lecavelier des Etangs et al., 
[20081 Hu itson et al.. l20T2t and super-Earth GJ 1214b (Berta et 
al., 120121 ) suggest that haze must be invoked to explain the ap- 
pearance of their transmission spectra. 
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Table 1. Total refracted angles c re f r for rays at 0.55 /urn entering the atmosphere with an impact radius rt,. For comparison, a* & is 
the refracted angle for an isothermal atmosphere of refractivity v^o^m an d scale height H atm (Baum & Code. [l953T l. 
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Fig. 1. Limb darkening function U(fi e ) at selected wavelengths. For a color version of the figure, see the electronic journal. 
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Fig. 2. Sketch of the Venus-Sun system at an instant near internal contact. The rays traced from x<j fan out as they traverse the 
Venus stratified atmosphere. 
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Fig. 3. The sketch demonstrates how a differential element of solid angle dQ. fans out when traced sunwards from the observer's site. 
The dashed line represents the trajectory passsing through the middle-point of the solid angle element. In a transparent atmosphere, 
the ratio of radiance and the squared refraction index is conserved on the ray trajectory, and the corresponding irradiance for the 
differential solid angle element is simply L(xo, So)so ■ n Xo c/Q(so)~ £(x , s Q )c/Q(So). For the calculation of Eq. Q, a mathematical 
transformation is applied to rewrite the integral in the the auxiliary variables £ and r\. 
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Fig. 4. Haze extinction coefficient at 60 km for cloud top at 70 km. At other altitudes, y(A; z)-y(A; 60 km)xp/p(6Q km). 
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Fig. 5. Scattering phase function for mode 1 (dashed) and mode 2 (solid) haze at selected wavelengths. A scattering angle of zero 
corresponds to the forward scattering direction. For a color version of the figure, see the electronic journal. 
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Fig. 6. Isolines of <D=- log 10 (Lfl(xo, So)/B Q ) at 0.55 jim during ingress and egress. The Venus solid core plus atmosphere up to 50 
km are represented to scale in dashed fill. The atmospheric ring from 50 to 100 km has been resized for a better appreciation of the 
halo's structure. The atmosphere above 100 km is not shown. The unrefracted view of the solar disk is shown to scale in yellow in 
the background. The directions for the Venus North Pole and the Celestial North are indicated. Also indicated is the azimuthal angle 
<p as defined in Figs. In the color code of the figure, blue indicates bright and brown indicates faint. For a color version of 

the figure, see the electronic journal. 
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Fig. 7. Limb-integrated lightcurves during ingress and egress at 0.43 (blue), 0.55 (green) and 0.8 (red) fim. Reference /! c i ou d-latitude 
profile. For a color version of the figure, see the electronic journal. 
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Fig. 8. Relative intensity against position angle from G-band observations (symbols) at the Dutch Open Telescope on La Palma 
(Spain) and modelling in this work at 0.43 /mi (solid) for the egress during the 2004 Venus transit. The intensities and position angle 
are defined by Tanga et al. d2012.l t. To reproduce the peak in brightness away from the South Pole, we needed to shift the reference 
^cloud-latitude profile by about 20°. For a color version of the figure, see the electronic journal. 
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Fig. 9. Limb-integrated lightcurves at 0.55 fi. Each pair of curves corresponds to the lightcurves obtained from perturbing the 
reference /i c i 0U( j-latitude profile by ± 1 km. 
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Fig. 10. Diffuse radiance L</(xo, So)/Bq graphed against the impact radius - Ry. From right to left, the curves correspond to 
wavelengths of 0.43, 0.55 and 0.8 ft, respectively. For a color version of the figure, see the electronic journal. 
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Fig. 11. Transit depths for the Venus-Sun integrated signal as observed from Earth. The dashed curves represent the differential 
transit depth contributed by refraction. 
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Fig. 12. Same as Fig. (fTTI) with the observer at an infinite distance. 
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Fig. 13. Mid-transit transmission spectrum of Venus at a resolving power of 10 3 in two different atmospheric conditions. For a color 
version of the figure, see the electronic journal. 
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